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Verfahren zur automatischen Ermittlunq optischer Parameter eines 

Schichtstapels 



Die Erfindung bezieht sich auf ein Verfahren zur automatischen Ermittlung 
5 optischer Parameter eines Schichtstapels, wie Schichtdicken, 

Brechungsindizes oder Absorptionskoeffizienten, durch Vergleich eines von 
einem Ort des Schichtstapels aufgenommenen optischen Mefcspektrums mit 
einem anhand vorgegebener optischer Parameterwerte berechneten 
Analysespektrum und Optimierung des berechneten Analysespektrums auf 
10 das Mefcspektrum hin. Die Erfindung bezieht sich weiterhin auf ein 

Computerprogramm (-produkt) zur Durchfuhrung eines solchen Verfahrens. 

Derartige Verfahren spielen insbesondere bei der Messung der Schichtdicke 
von dunnen Schichten und weiterer optischer Parameter, wie Brechungsindex 
15 und Extinktionsfaktor, von Ein- und Mehrschichtsystemen, wie sie 
beispielsweise strukturierte Wafer darstellen, eine wichtige Rolle. 

In vorliegender Beschreibung umfasst der Begriff "Schichtstapel" sowohl den 
Schichtstapel im engeren Sinne (Abfolge einzelner Schichten, wie Si02, 
20 Si3N4, Resist-Filme, etc. auf einem Substrat, wie Silizium oder Aluminium) als 
auch die Kombination aus Schichtstapel (Layer) und Substrat. 

Eine optische Messeinrichtung zur Messung der genannten Eigenschaften an 
Ein- und Mehrfachschichtsystemen in einem Schichtdickenbereich von etwa 1 
25 nm bis zu etwa 50 pm ist aus der DE 100 21 379 A1 bekannt. Vorgesehen ist 
dort eine Beleuchtungseinrichtung, beispielsweise eine Halogen- und eine 
Deuteriumlampe, um einen Messlichtstrahl eines ausreichend breiten 




Wellenlangenbereichs, beispielsweise zwischen 190 nm und 800 nm, zu 
erzeugen. Mittels eines Strahlteilers wird der Messlichtstrahl in einen 
Objektlichtstrahl und einen Referenzlichtstrahl aufgeteilt. Der Messlichtstrahl 
wird mittels eines Messobjektivs in annahernd senkrechtem Einfall auf den 
Messort einer Probe gelenkt und der von der Probe reflektierte Strahl wird 
zusammen mit dem Referenzlichtstrahl einer Auswerteeinrichtung zugefuhrt. 
Eine geeignete Auswerteeinrichtung ist hierbei ein Spiegel-Gitter- 
Spektograph, der die Wellenlangen des auftreffenden Lichtes raumlich 
getrennt auf einen CCD-Detektor abbildet. Dieser ist im gesamten 
Wellenlangenbereich empfindlich und erlaubt ein schnelles Auslesen der 
Messspektren. In der genannten Schrift wird der reflektierte Objektlichtstrahl 
sowie der Referenzlichtstrahl uber Lichtleiter der Auswerteeinheit zugefuhrt. 
Die beschriebene Messanordnung kann zusatzliche eine einkoppelbare 
Einrichtung zur visuellen Darstellung und Uberwachung enthalten. 

Mit einer Messanordnung gemafJ DE 100 21 379 A1 lassen sich zur 
Ermittlung der optischen Schichteigenschaften die aufgrund Interferenzen 
entstehenden Intensitatswerte im Spektrum des von der Probe reflektierten 
Objektlichtstrahls detektieren und auswerten. Wegen der Vieldeutigkeit (die 
Intensitatswerte berechnen sich je nach Schichtfolge aus mehreren Termen, 
die eine Funktion zum Sinus der Phase aus dem Produkt von jeweiliger 
Schichtdicke und (spektral abhangigen) Brechungsindex, sowie den 
Brechungs- und Absorptionsindizes selbst sind) konnen, von Sonderfallen 
abgesehen, aus der Kurvenform die optischen Parameter analytisch nicht 
zuruckgerechnet werden. In der Regel mussen rechenintensive Fitverfahren 
eingesetzt werden. 

Aus dem Stand der Technik sind verschiedene Verfahren zur Auswertung des 
Spektrums des reflektierten Objektlichtstrahls bekannt. Beispielsweise lasst 
sich gemali der europaischen Patentschrift EP 0 644 399 B1 die Schichtdicke 
d einer diinnen Einfachschicht aus der Anzahl m der Extremwerte (Maxima 
und Minima) im Spektrum des reflektierten Objektlichtstrahls im beobachteten 
Wellenlangenbereich von A.1 bis X2 aus der bekannten Formel 
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bestimmen, wobei n1 und n2 jeweils der Brechungsindex der dunnen Schicht 
bei der Wellenlange A.1 bzw. X2 ist. 

5 

Bei Mehrfachschichtsystemen erhalt man jedoch ein Spektrum, in dem die 
Interferenzspektren der einzelnen Schichten sowie der Schichten 
untereinander Ciberlagert sind, so dass die Gleichung (1) nicht mehr 
unmittelbar anwendbar ist. In einem solchen Fall konnen globale und lokale 

10 Optimierverfahren eingesetzt werden, die von theoretischen Modellen mit 
vorgegebenen Schichtdickenbereichen ausgehen und diese hinsichtlich des 
ermittelten Spektrums optimieren. Das Verfahren gemafl der genannten 
Patentschrift geht von einem mdglichen Schichtdickenbereich aus, der von der 
Gesamtzahl der Extrema, der Wellenlange des untersten und des obersten 

15 Extremums sowie einem uber den Wellenlangenbereich gemittelten 

Brechungsindex einer Schicht abhangt Durch Veranderung der Schichtdicke 
im jeweiligen Schichtdickenbereich mit vorbestimmten Schrittweiten fur jede 
einzelne Schicht wird die Schichtdickenkombination bestimmt, deren 
berechnete spektrale Reflexion die geringste Abweichung zur gemessenen 

20 aufweist. 

Das Verfahren der EP 0 644 399 B1 stellt kein allgemeines Verfahren mit 
Variationsmoglichkeiten von Brechungs- und Absorptionsindex dar, da diese 
optischen Eigenschaften einer jeden Schicht sowie die Anzahl der Schichten 
25 bekannt sein mussen. Die Schichtdickenbereiche besitzen immer 0 als untere 
Grenze, ausgewertet werden nur die Extrema-Lagen. 

Die US-4,984,894 misst die Dicke der obersten Schicht eines 
Mehrschichtsystems unter der Voraussetzung, dass von der darunter 
30 liegenden zweiten Schicht kein Licht reflektiert wird. 



Das genannte Verfahren ist auf die oberste Schicht einer bestimmten 
Schichtenfolge sowie auf bestimmte Schichtenparameter eingeschrankt und 
lasst nur naherungsweise Ergebnisse zu. 

In der US-5,440,141 werden die Schichtdicken eines 3-fach-Schichtsystems 
bekannter Zusammensetzung bestimmt, indem fur die oberste Schicht die 
bereits behandelte Extrema-Methode und fur die beiden folgenden Schichten 
eine Fburiertransformations-Methode zusammen mit Optimierungsverfahren 
fur die erhaltenen Schichtdicken eingesetzt werden. Bei der Fourier-Methode 
wird das in Abhangigkeit der Wellenlange gemessene Reflexionsspektrum in 
ein von der Wellenzahl abhangiges Spektrum konvertiert und anschlielSend 
Fourier-transformiert. Der Absolutbetrag des Fourier-transformierten 
Spektrums zeigt im Falle einer 2-fach-Schicht drei Peaks, einen fur jede 
Schicht und einen Summenpeak. Diese Peaks genugen der Summenrelation, 
so dass unpassende Peaks ausgeschlossen werden konnen. Einem Peak im 
Fourier-transformierten Spektrum lasst sich ein Wert der optischen Dicke (nd) 
zuordnen, wenn die optischen Dicken in Relation zum gemessenen 
Spektralbereich genugend dick (mindestens eine Periode im Spektralbereich) 
sind. 

Das geschilderte Verfahren der US-5, 440,141 ist auf bestimmte 
Schichtenkombinationen bekannter Zusammensetzung eingeschrankt und 
lasst sich nicht fur allgemeine Messungen heranziehen. 

Schliefilich ist aus der US-5,864,633 ein Verfahren zur optischen Inspektion 
eines Filmstacks (Dunnschichtstapel) bekannt, bei dem optische Daten und 
dazu korrespondierende theoretische Daten verglichen und die theoretischen 
Daten mittels genetischer Algorithmen angepasst werden. Jedes theoretische 
Modell reprasentiert hierbei einen sogenannten Genotyp (Satz von 
Dunnschichtparametern), der eine geordnete Liste von Genen (verschiedene 
Schichtparameter wie Dicke, Brechungsindex, Extinktionskoeffizient) darstellt. 
Ein Genotyp enthalt somit die verschiedenen Schichtparameter samtlicher 
Schichten. Zunachst wird eine Anzahl von Genotypen definiert und fur jeden 



Genotyp ein Fit-Level bestimmt, der sich aus dem Vergleich der berechneten 
theoretischen Daten mit Messdaten der optischen Inspektion ergeben. Je 
nach Fit-Level werden die Genotypen einer genetischen Operation (Kopieren, 
Kreuzung, Mutation) unterzogen. Aus dem vorhandenen Satz von Genotypen 
kann auf diese Weise ein neuer Satz von Genotypen (neue Generation) 
erzeugt werden. Verbessert sich der Fit-Level des besten Genotyps uber eine 
Anzahl von Generationen nicht mehr wesentlich, so wird das Verfahren 
abgebrochen. 

Aufgrund der hohen Anzahl von Rechenoperationen und der daraus 
resultierenden Rechenzeit ist dieses Verfahren fur den industriellen Einsatz 
der Uberprufung und Vermessung von Schichtensystemen nicht geeignet. 

Der vorliegenden Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren 
zur automatischen Ermittlung optischer Materialeigenschaften eines 
Schichtstapels anzugeben, das ohne Einschrankungen hinsichtlich Anzahl, 
Beschaffenheit oder Dicke der Schichten mit moglichst wenig 
Rechenoperationen und somit in kurzer Zeit Ergebnisse liefert, die den 
Einsatz dieses Verfahrens insbesondere in der kontinuierlichen 
Produktionslinie, beispielsweise bei der Wafer-Fertigung, erlauben. 

Diese Aufgabe wird durch ein Verfahren der eingangs genannten Art gelost, 
bei dem das aufgenommene Meftspektrum anhand von 
Kurvenformparametern, die das MeSspektrum charakterisieren und aus 
diesem ermittelt werden, klassifiziert wird und diese Kurvenformparameter mit 
entsprechenden fur bekannte Schichtstapel berechnete Kurvenformparameter 
von Spektren verglichen werden, urn Werte oder Wertebereiche fur die zu 
bestimmenden optischen Parameter zu ermitteln, anhand derer das oder die 
Analysespektren zum Vergleich mit dem MeBspektrum berechnet werden. 
Aus der Liste der berechneten Kurvenformparameter werden relevante Werte 
oder Wertebereiche durch den Vergleich mit den gemessenen heraus gefiltert. 
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Die erfindungsgemafte Klassifizierung des aufgenommen MeUspektrums 
anhand charakteristischer Kurvenform para meter und der anschliefiende 
Vergleich mit entsprechenden fur bekannte Schichtstapel berechneten 
Kurvenformparametern fuhrt unmittelbar zu einem ersten Ergebnis fur die zu 
5 bestimmenden optischen Parameter. Hieraus wird das Analysespektrum 
berechnet und mit dem gemessenen Meftspektrum verglichen. Je nach 
Qualitat der Ubereinstimmung schlielien sich hieran weitere Fit-Verfahren, wie 
weiter unten erlautert wird, an. Der genannte Vergleich der 
Kurvenformparameter der klassifizierten optischen Spektren kann auch 
10 Wertebereiche fur die zu bestimmenden optischen Parameter ergeben, die 
den nachfolgenden Fit-Verfahren zugrunde gelegt werden. 



Entscheidender Vorteil des erfindungsgemafJen Verfahrens ist die 
Einschrankung des mogtichen Wertebereichs fur die zu bestimmenden 

15 optischen Parameter eines Schichtstapels mittels Vergleich von 

Spektrenparametern (Kurvenformparameter), wobei dieser Vergleich durch 
Verwendung von zuvor berechneten und vorsortierten Tabellen in 
vergleichsweise kurzer Zeit bewerkstelligt werden kann. Fur nachfolgende Fit- 
Verfahren stehen folglich im Vergleich zu bisherigen Methoden wesentlich 

20 reduzierte Wertebereiche fur die zu bestimmenden optischen Parameter zur 
Verfugung, so dass diese Fit-Verfahren in wesentlich kurzerer Zeit zu 
bewerkstelligen sind. 

i 

Im folgenden sei die Erfindung einer herkommlichen Naherungsmethode, wie 
25 sie in der bereits erwahnten EP 0 644 399 B1 behandelt ist, am Beispiel einer 
Dreifachschicht gegenubergestellt. Zum Einsatz kommt die Untersuchung des 
Reflexionsspektrums im Bereich von 400 bis 800 nm. Die Gesamtdicke sei 
derart, dass mehrere Extrema entstehen: 



30 



a) Zunachst werden alle 401 Spektrenkanale ausgewertet und die Anzahl 
der Extrema bestimmt; 



b) Die Schichtdickenabschatzung gemafi genanntem Patent ergebe 
obere Grenzwerte von 700nm, 500nm und 400nm; 

c) Fur den Grob-Fit werde die Schichtdicke im Bereich von Null bis zum 
jeweiligen oberen Grenzwert in Schrittweiten von 10nm variiert; 

d) Es ergeben sich 70 x 50 x 40= 140.000 Stutzstellen fur die 
Dickenberechnung, das heilit es mussen 140.000 Spektren berechnet 
und verglichen werden. Daran anschlieliend folgt der sogenannte 
Fein-Fit, indem in einem weiteren Iterationsverfahren das lokale 
Minimum genau bestimmt wird. Auch hierbei wird in jedem 
Iterationsschritt ein Theoriespektrum berechnet. 

Beim geschilderten Beispiel sollten als optische Parameter nur die 
Schichtdicken bestimmt werden. Weitere Parameter wie Brechungsindex oder 
Absorptionskoeffizient gehen beim Grob- als auch beim Fein-Fit multiplikativ 
ein, so dass die Anzahl der Stutzstellen schnell mehrere Millionen betragen 
kann. Die typischen Auswertezeiten liegen bei dem genannten Beispiel 
oberhalb derjenigen, die das Verfahren fur den kontinuierlichen industriellen 
Einsatz tauglich machen. 

Die beispielhaft beschriebene Methode des genannten Patents hat weiterhin 
den Nachteil, dass als untere Grenze immer der Wert Null fur die Schichtdicke 
angenommen werden muss, wenn keine weiteren Beschrankungen 
vorgegeben sind. Wird der Parameterraum zu stark eingeschrankt, urn die 
Analysezeit zu verkurzen, so kann es zu Fehlauswertungen kommen. Wird 
der Parameterraum zu grob nach lokalen Minima abgesucht, besteht ein 
grolies Restrisiko hinsichtlich der Fehlinterpretation der Daten und einer 
Auswertung, die zu einem falschen Ergebnis fuhrt Des weiteren kann durch 
die Interferenzeffekte die Bestimmung der Anzahl der Extrema falsch sein, ein 
Fehler, der sich bei der Bestimmung der oberen Grenzwerte fur die 
Schichtdicken auswirkt und entsprechend fortpflanzt. 
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Eine Einschrankung des Parameterraumes ist aus Zeitgrunden bei der 
Auswertung dringend erforderlich, insbesondere wenn die Anzahl der 
Schichten ansteigt. 

Gemafc vorliegender Erfindung wird das aufgenommene MelSspektrum mittels 
charakteristischer Kurvenformparameter klassifiziert, wobei in der Regel in der 
Grolienordnung 5 bis 15 solcher Parameter ausreichend sind. Die 
Kurvenformparameter des aufgenommenen MeUspektrums werden nun mit 
den tabellierten Kurvenformparametern bekannter Spektren verglichen, wobei 
man als Ergebnis fur jeden zu bestimmenden optischen Parameter einzelne 
Werte Oder einen Wertebereich erhalt. Folglich werden bei der Erfindung 
zunachst nicht Spektren im Umfang von 400 bis 600 Werten miteinander 
verglichen, sondern Tabelleneintrage (von ca. 10 Werten), womit eine 
erhebliche Reduktion von Rechenkapazitat und -zeit verbunden ist. 
Ausschlaggebend bei der Zeitersparnis ist, dass die Berechnung eines 
Spektrums uber eine komplexe Formel eine wesentlich grossere Zeit (> Faktor 
100000) in Anspruch nimmt als der Vergleich mit den Tabelleneintragen. 

Die Klassifizierung des MelJspektrums erfolgt anhand eines oder mehrerer der 
folgenden charakteristischen Kurvenformparameter: das lokale Rauschen 
des Spektrums, der Mittelwert, die Standardabweichung des Mittelwerts, die 
Anzahl und Lage der Extrema, eine Klassifizierung der Extrema, 
beispielsweise nach spektraler Lage, den Intensitatswerten oder den relativen 
Abstanden zueinander, Parameter der einhullenden Kurven der Minima und 
Maxima, der gemittelte Kurvenverlauf, Schwebungen sowie mogliche weitere 
Parameter, wie die Anzahl der Peaks im Fourier-transformierten Spektrum. 

Eine Einschrankung bzw. Filterung der Wertebereiche fur die zu 
bestimmenden optischen Parameter erfolgt vorzugsweise durch Vergleich der 
Spektrenparameter mit vorgefertigten Parameterlisten (Tabellen) und je nach 
Schichtstapel zusatzlich noch durch eine Extremamethode und/oder durch 
eine Fourier-Transformations-Methode. Beispiele solcher Methoden sind aus 
dem Stand der Technik - wie eingangs erwahnt - bekannt. 
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Die Ermittlung der optischen Parameter des untersuchten Schichtstapels kann 
auf der Grundlage des eingeschrankten Parameterraumes vorteilhaft 
anschlieRend mit bekannten Grob- und Fein-Fit- Verfahren, beispielsweise 
5 mittels Raster, Intervallmethode- und/oder Powell-Methode, erfolgen. Die 
Ubereinstimmung von Mefi- und Analysespektrum wird daraufhin bewertet 
und der "best fit" ausgewahlt. 

Sollte die beschriebene Methode nicht zu plausiblen Ergebnissen fuhren, kann 
10 gegebenenfalls der eingeschrankte Parameterraum erweitert und das 
Verfahren von Neuem durchlaufen werden. 

Haufig ist der Aufbau des Schichtstapels, das heillt die Abfolge der 
Zusammensetzung der einzelnen Schichten bekannt. Andernfalls kann in 

1 5 einer Ausgestaltung des erfindungsgemaBen Verfahrens, fur die gesondert 
Schutz beansprucht wird, in einem ersten Schritt eine automatische 
Bestimmung der Abfolge der Zusammensetzung des Schichtstapels erfolgen, 
indem wiederum ein MeUspektrum aufgenommen und anhand 
charakteristischer Kurvenformparameter klassifiziert wird und durch Vergleich 

20 mit entsprechenden Kurvenformparametern von zu Schichtstapeln bekannter 
Zusammensetzung gehorenden Spektren eine oder mehrere mdgliche 
Abfolgen der Schichtstapelzusammensetzung bestimmt werden. 

Weiterhin kann auch in diesem Fall anhand der Ergebnisse der 
25 Schichtstapelzusammensetzung ein Analysespektrum berechnet werden, das 
mittels Fit-Verfahren auf das Meftspektrum hin optimiert wird. Gleichzeitig 
konnen neben einer moglichen Abfolge der Schichtstapelzusammensetzung 
auch die Schichtdicken-Bereiche, Brechungsindex-Bereiche sowie weitere 
Bereiche fur die in Frage kommenden optischen Parameter bestimmt werden. 
30 In diesem Fall muss in einem wesentlich groReren Parameterraum gesucht 
werden, so dass es vorteilhaft ist, diese Vorabbestimmung der 
Schichtstapelzusammensetzung und seiner optischen Parameter im 
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Hintergrund, beispielsweise simultan zum Einlernen der Tischpositionen, 
vorzunehmen. 

Haufig kann der zu durchsuchende Spektrenparameterraum eingeschrankt 
werden, indem der Kunde die bei ihm zum Einsatz kommenden moglichen 
Layer-Substrat-Kombinationen vorgibt. Aus den vorhandenen Moglichkeiten 
sucht das erfindungsgemaUe Verfahren dann vorab die wahrscheinlichsten 
Kombinationen (sowie die zugehorigen optischen Parameterbereiche) heraus. 

Es ist vorteilhaft, die in dieser Vorabbestimmung gefundenen Ergebnisse dem 
Kunden anzuzeigen und ihm die Moglichkeit zu geben, das Ergebnis zu 
ubernehmen oder zu korrigieren. 

Die erfindungsgemafJe Ermittlung optischer Parameter eines Schichtstapels 
mit der etwaigen Ermittlung der Abfolge der chemischen Zusammensetzung 
des Schichtstapels wird in vorteilhafterweise mittels eines 
Computerprogramms durchgefuhrt, das auf einer geeigneten Rechnereinheit 
ausgefuhrt wird. Ermittelte Daten (Wertebereiche fur die optischen Parameter, 
Schichtzusammensetzung) konnen in gangiger Weise auf einem Monitor 
angezeigt werden. Weiterhin kann die Moglichkeit einer Beeinflussung der 
angezeigten Daten dem Kunden eingeraumt werden. Das Computerprogramm 
kann auf geeigneten Datentragern, wie EEPROMs, Flash-Memories, aber 
auch CD-ROMs, Disketten oder Festplattenlaufwerken gespeichert sein. Auch 
eine Ubertragung des Computerprogramms uber ein Kommunikationsmedium 
(wie das Internet) zum Kunden (Anwender) ist moglich. 

Im folgenden soil anhand der beigefugten Figuren ein Ausfuhrungsbeispiel die 
Erfindung und deren Vorteile naher erlautern. 
Es zeigt 

Figur 1 Zwei MeBspektren einer Zweifachschicht (Fig. 1a) und einer 
Einfachschicht (Fig. 1b) auf einem Substrat, 
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Figur2 den Kurvenformparameter der Anzahl der Extrema aufgetragen uber 
die optische Dicke einer berechneten Parameterliste fur die genannte 
Zweifach-Schicht, 

5 Figur 3 den Mittelwert des Spektrums als Kurvenformparameter aufgetragen 
uber die optische Dicke einer fur die genannte Zweifach-Schicht berechneten 
Parameterliste, 

Figur 4 den Wellenlangenwert des Maximums, der dem langwelligen Ende 
10 der Meftspekrums am nachsten liegt, des berechneten Spektrums als 
Kuvenform para meter aufgetragen uber die optische Dicke der fur die 
genannte Zweifach-Schicht berechneten Parameterliste, 

Figur 5 der Maximalwert als Kurvenformparameter aufgetragen uber die 
15 optische Dicke der fur die genannte Zweifach-Schicht berechneten 
Parameterliste, 

Figur 6 theoretische Spektren ahnlichen Aussehens 

20 Die Erfindung soil im Folgenden anhand des einfachen Beispiels einer 
Zweifach-Schicht auf einem Substrat erlautert werden, sie ist jedoch 
keinesfalls auf diesen Spezialfall beschrankt. Das Beispiel bedient sich einer 
Si3N4-Si02-Si-Kombination (Si als Substrat). Fur das nachfolgend 
geschilderte erfindungsgemafce Verfahren ist folglich in diesem Fall die 

25 Abfolge der Zusammensetzung des Schichtstapels bekannt. 

Die unten stehende Tabelle stellt ein Beispiel einer vorgefertigten berechneten 
Parameterliste fur die genannte Zweifach-Schicht dar, wobei zu vorgegebenen 
Dickenwerten D1 (Dicke der Schicht Si3N4) und D2 (Dicke der Schicht Si02) 
30 und der daraus resultierenden optischen Gesamtdicke die darauffolgenden 
Kurvenformparameter aus den zugehorigen berechneten Analysespektren 
abgeleitet wurden: 
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Vorgegebene Dicke der ersten Schicht 



"D2" 
5 "Opt.Thick." 



Vorgegebene Dicke der zweiten Schicht 

Optische Dicke (Optical Thickness) berechnet aus der 
Summe der Produkte von mittlerem Brechungsindex 
und Dicke im Wellenlangenbereich von 200 nm bis 800 
nm, also Optische Dicke = <n1(A.)> D1 + <n2(A.)> D2 



10 "NoE" 



Anzahl der Extrema (Number of Extrema) 



15 



"Mean" 



"Sigma" 



Mittelwert des berechneten Spektrums 

Standardabweichung zum Mittelwert des berechneten 
Spektrums 



20 



"Min" 



"Max" 



"WL-MaxEx" 



Intensitatswert des Minimums im Wellenlangenbereich 

Intensitatswert des Maximums im Wellenlangenbereich 

Wellenlange, bei der das letzte Maximum auftritt, 
beginnend bei der kleinsten Wellenlange <in nm> 



25 



"MDEx" 



Mittlerer Abstand der Extrema (= "MeanDistExtrema") 
bei mehr als einem Extremum <in nm> 



Weitere sinnvolle Kurvenformparameter waren, insbesondere bei dickeren 
Schichten, die sich aus einer Fast-Fourier-Transformation ergebenden Werte, 
wie Lage der Einzelpeaks und des Summenpeaks. Ebenso suchen konnte 
30 man nach auftretenden Schwebungen Oder nach Lage und Intensitat der 
auftretenden Extrema 
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Nachfolgend wiedergegeben ist die der Auswertesoftware als Look-Up- 
Tabelle vorliegende Kurvenformparameterliste. 

Tabelle 1 



Si3N4Si02Si 



Input 




Output 
















D1 


D2 


Opt. 
Thick 


NoE 


Mean 


Sigma 


Min 


Max 


WI-MaxEX 


MDEx 


0 


0 


0 


0 


1.000 


0.000 


1.000 


1.000 


0 


0 


0 


20 


30 


0 


0.909 


0.071 


0.663 


0.971 


0 


0 


20 


0 


43 


1 


0.752 


0.192 


0.264 


0.933 


245 


0 


0 


40 


60 


1 


0.753 


0.091 


0.559 


0.892 
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Die in Figur 1 dargestellten Melispektren konnen beispielsweise mit einer 
optischen Messeinrichtung aufgenommen werden, wie sie aus der eingangs 
behandelten DE 100 21 379 A1 bekannt ist. Zu Einzelheiten der Messung wird 
voll umfanglich auf diese Schrift verwiesen. Aus dem aufgenommen 
MefSspektrum werden erfindungsgemali die in diesem Beispiel genannten 
charakteristischen Kurvenformparameter abgeleitet und die Ergebnisse mit 
den Werten der oben stehenden Tabelle verglichen. Als Ergebnis erhalt man 
nun eine oder mehrere optische Dicken und somit Schichtdicken- 
Kombinationen, fur die eine besonders gute Ubereinstimmung der aus dem 
Mefcspektrum abgeleiteten Kurvenformparameter mit den berechneten 
Parametern der Liste gegeben ist. Fur diese Dickenkombination werden nun 
zugehorige Analysespektren berechnet, die mit einem aufgenommenen, wie in 
Figur 1a dargestellten, Mefcspektrum verglichen werden. Da in der Regel nicht 
davon auszugehen ist, dass die mit dem erfindungsgemafien Verfahren 
gefundene Dickenkombination bereits der vorhandenen entspricht, schlie&en 
sich zur Ermittlung der genauen Schichtdicken vorteilhafterweise bekannte 
Grob- und Fein-Fit-Verfahren, wie Raster-, Intervall- und Powell-Methode an. 
In diesem Fall dient das erfindungsgemafce Verfahren zur Einschrankung des 
Parameterraums, so dass die anschliefienden Fit-Verfahren wesentlich 
schneller zum Ziel fuhren. 

Es ist vorteilhaft, neben dem erfindungsgemaden Verfahren zur 
Einschrankung des Parameterraumes weitere bekannte Methoden 
hinzuzuziehen, insbesondere um beispielsweise aufgefundene Schichtdicken- 
Kombinationen (D1, D2) als unplausibel auszuschlie&en. Hierzu konnen 
insbesondere die bereits erwahnte Extrema-Methode und die Fourier- 
Transformations-Methode verwendet werden. 
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Die Figuren 2 bis 5 zeigen, wie bestimmten charakteristischen 
Kurvenformparametern der aufgenommenen Messkurve vorgegebene Werte 
von optischen Parametern (in diesem Fall Schichtdicken-Kombinationen) 
zugeordnet werden konnen. In Figur 2 ist der Zusammenhang angenahert 
linear, d.h. die Anzahl der Extrema (Number of Extrema) nimmt proportional 
zur optischen Dicke (Opt. Thickness) zu. Der Parameter Mittelwert (Mean) 
(Figur 3) andert sich in Form einen gedampften Schwingung mit der optischen 
Dicke, wobei die Schwankungsbreite mit zunehmender optischer Dicke zwar 
abnimmt, der Mittelwert sich jedoch auch immer mehr einer Konstanten 
nahert. Dies spiegelt naturlich auch die spektrale Auflosung der MeGapparatur 
und damit das Abtasttheorem wider. Der in Figur 4 uber die optische Dicke 
aufgetragene Parameter WLMaxEx beschreibt die Lage des langwelligsten 
Maximums. Diese Werte sind naturlich auch durch den Wellenlangenbereich 
der MeUapparatur (hier 200nm und 800nm) begrenzt. Kurven die kein 
eindeutiges Maximum zeigen (die Randwellenlangen werden ausgeschlossen, 
Extrema mussen einen vordefinierten Schwellenwert ubersteigen), erhalten 
als Parameterwert Null zugeordnet. Von der optischen Dicke Null ausgehend 
wachst dieser Wert an, bis (sozusagen) das Extremum aus dem Mefcbereich 
herausgewandert ist. Figur 5 zeigt, dass der Parameter Maximum (Maximum 
Value)(lntensitatswert) in dem angezeigten Bereich der optischen Dicke 
annaherend von Wert zu Wert oszilliert. 

Insgesamt folgt, dass die Zuordnung einer optischen Dicke uber einen 
einzelnen aus der Melikurve gewonnenen Wert vieldeutig ist. Daher mussen 
mehrere dieser Werte herangezogen werden. Die Schwankungsbreiten in den 
einzelnen Kurven zeigen an, wie verschieden stark die Wertebereiche der 
Parameter einzuengen sind. Die Moglichkeit der Einengung und damit der 
Filterung ist durch die horizontalen Linien in den Figuren als Beispiel einer 
mogliche Auswertevariante verdeutlicht. 

Allgemein gilt: Aus den zugeordneten Werten lassen sich die 
wahrscheinlichsten auswahlen (auch mit Hilfe weiterer bekannter Methoden) 
und zur Berechnung eines Analysespektrums verwenden. 
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Vereinfachtes Ausfuhrungsbeispiel Filter: 

Aus einem Spektrum (175nm Si3N4 auf 190nm Si02), das nicht dem in Fig. 
1a dargestellten entspricht, ergeben sich die in Tabelle 2, Tabellenspalte 1 
„Gesuchter Wert" stehenden Such-Werte. 

Selektiert man aus den ursprunglich in Tabelle 1 angegebenen 256 
Listeneintrage nacheinander die in der Tabelle 2 angegebenen Filterbereiche 
(dies entspricht den horizontalen Linien in den Figuren 2-5), so reduziert sich 
sukzessive die Anzahl der Listeneintrage von zunachst 63 auf 4. 

Die diesen Listeneintragen zugeordneten Spektren sind in Fig. 6 zusammen 
mit dem gesuchten Spektrum (175 - 190) dargestellt. 

Die beste Ubereinstimmung ergibt sich fur die benachbarten Kurven mit den 
Schichtdicken (180-180) und (160-200). 



Tabelle 2: Reduktion durch Filtern. Ausgangswert 256 Listeneintrage 



Gesuchter 


Filterstufe/ 


Filter-bereich 


Anzahl der 


Wert 


Filtername 




Listeneintrage 


11 


1 /NoE 


10-12 


63 


628 


2 / WL-Max 


613-643 


12 


1.4 


3 /Max 


1.25-1.55 


9 


0.78 


4 / Mean 


0.74-0.82 


4 



Ein Grobfit in den angegebenen Dickenbereichen (z. B. Schrittweite jeweils 
±20nm) wird zu einem Ergebnis mit guter Ubereinstimmung der Kurvenform 
fur den Tabelleneintrag mit den Dicken D1 = D2 = 180nm fuhren. Ein 
anschliessender Feinfit mittes Raster-, Intervall- oder Powell-Methode fuhrt zu 
einem Ergebnis der gewunschten Genauigkeit (z.B. 0.1 nm). 



A 3041 DE 
11.07.2002 
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Aus Grunden der Einfachheit beschrankt sich das genannte Beispiel auf die 
Bestimmung lediglich der Schichtdicken einer Zweifach-Schicht. Dem 
Fachmann ist verstandlich, wie das Beispiel auf die Bestimmung weiterer 
optischer Parameter, wie Brechungsindex n oder Extinktionskoeffizient k, 
5 ausgedehnt werden kann. 



v. 



Insbesondere ist mit dem genannten erfindungsgemaften Verfahren auch eine 
Vorauswahl der in Frage kommenden Schichtarten (chemische 
Zusammensetzung) moglich, wobei hier aus einem entsprechend grolieren 



10 Parameterraum (Parameterlisten fur verschiedene Zusammensetzungen von 
Ein- oder Mehrfach-Schichten) gewahlt werden muss. Eine Vorab- 
Beschrankung ist jedoch meist moglich, da dem Kunden (Anwender) meist 
bekannt ist, urn welche moglichen Kombinationen es sich handeln kann. 
Beispielsweise kann die Bestimmung der vorhandenen Kombination (also der 
15 Abfolge der Schichtzusammensetzungen) wahrend des Einlernens der 

Tischpositionen zur nachfolgenden Messung im Hintergrund erfolgen. Vor 
dem Beginn der eigentlichen Messung der optischen Parameter wird dann 
dem Kunden (Anwender) die gefundene Kombination angeboten, die er 
ubernehmen oder korrigieren kann. 
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Patentanspruche 



Verfahren zur automatischen Ermittlung optischer Parameter eines 
Schichtstapels, wie Schichtdicken, Brechungsindizes oder 
Absorptionskoeffizienten, durch Vergleich eines von einem Ort des 
Schichtstapels aufgenommenen optischen Meftspektrums mit einem 
anhand vorgegebener optischer Parameterwerte berechneten 
Analysespektrum und Optimierung des berechneten Analysespektrums 
auf das Melispektrum hin, 
dadurch gekennzeichnet, 
dass das aufgenommene MelJspektrum anhand von 
Kurvenformparametern, die das MeBspektrum charakterisieren und aus 
diesem ermittelt werden, klassifiziert wird und diese Kurvenformparameter 
mit entsprechenden fur bekannte Schichtstapel berechnete 
Kurvenformparametern von Spektren verglichen werden, urn Werte oder 
Wertebereiche fur die zu bestimmenden optischen Parameter zu 
ermitteln, anhand derer das oder die Analysespektren zum Vergleich mit 
dem Meftspektrum berechnet werden. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das 
aufgenommene Melispektrum anhand eines oder mehrerer der 
folgenden Kurvenformparameter klassifiziert wird: das lokale Rauschen 
des Spektrums, der Mittelwert des Spektrums, die Standardabweichung 
des Mittelwerts, die Anzahl und Lage der Extrema, eine Klassifizierung 
der Extrema, beispielsweise nach spektraler Lage, den Intensitatswerten 
oder den relativen Abstanden zueinander, Merkmale der einhullenden 
Kurven der Minima und Maxima, ein gemittelter Kurvenverlauf, 
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Schwebungen sowie Parameter aus der Fourier-Transformierten des 
aufgenommenen Mefcspektrums, wie Anzahl, Lage und Werte der dort 
vorhandenen Extrema. 

Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass zur 
Einschrankung der Wertebereiche fur die zu bestimmenden optischen 
Parameter je nach Art des Schichtstapels zusatzlich eine Auswertung des 
aufgenommenen Mefispektrums nach einer Extremamethode und/oder 
einer Fourier-Transformations-Methode erfolgt. 

Verfahren nach einem der Anspruche 1 bis 3, dadurch gekennzeichnet, 
dass die Optimierung des berechneten Analysespektrums auf das 
Melispektrum hin mittels bekannter Grob- und/oder Feinfitverfahren 
vorgenommen wird. 

Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, 
dass die fur die Optimierung des berechneten Analysespektrums 
ermittelten Werte fur die zu bestimmenden optischen Parameter 
gegebenenfalls korrigiert werden. 

Verfahren zur automatischen Ermittlung der Abfolge der 
Zusammensetzung eines Schichtstapels, insbesondere gemafc einem 
Verfahren nach einem der vorhergehenden Anspruche, 
dadurch gekennzeichnet, 

dass ein von einem Ort des Schichtstapels aufgenommenes optisches 
MeUspektrum anhand von Kurvenformparametern, die das MefJspektrum 
charakterisieren und aus diesem ermittelt werden, klassifiziert wird und 
durch Vergleich mit entsprechenden Kurvenformparametern von zu 
bekannten Schichtstapeln gehorenden klassifizierten Spektren eine oder 
mehrere mogliche Abfolgen der Zusammensetzung des Schichtstapels 
bestimmt werden. 
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7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, dass gleichzeitig 
zur Bestimmung der Abfolge der Zusammensetzung des Schichtstapels 
aus dem Vergleich der Kurvenformparameter der klassifizierten Spektren 
Wertebereiche fur die zu bestimmenden weiteren optischen Parameter 
ermittelt werden. 

8. Verfahren nach Anspruch 6 oder 7, dadurch gekennzeichnet, dass 
aufgrund der bestimmten Abfolge der Zusammensetzung des 
Schichtstapels sowie etwaiger weiterer optischer Parameterwerte 
Analysespektren berechnet werden, die auf das aufgenommene hin 
optimiert werden. 

9. Verfahren nach einem der Anspruche 6 bis 8, dadurch gekennzeichnet, 
dass die ermittelte Abfolge der Zusammensetzung des Schichtstapels 
sowie etwaige weitere ermittelte optische Parameter einer Uberprufung 
unterzogen werden, bevor die automatische Ermittlung optischer 
Parameter des Schichtstapels gema(J Anspruch 1 erfolgt. 

10. Computerprogramm mit Programmcode-Mitteln, urn alle Schritte eines 
Verfahrens gemali einem der Anspruche 1 bis 9 auszufuhren, wenn das 
Computerprogramm auf einem Computer oder einer entsprechenden 
Rechnereinheit ausgefuhrt wird. 

11. Computerprogramm-Produkt mit Programmcode-Mitteln, die auf einem 
computerlesbaren Datentrager gespeichert sind, urn ein Verfahren nach 
einem der Anspruche 1 bis 9 auszufuhren, wenn das Computerprogramm 
auf einem Computer oder einer entsprechenden Rechnereinheit 
ausgefuhrt wird. 
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Zusammenfassung 



Die Erfindung betrifft ein Verfahren zur automatischen Ermittlung 
optischer Parameter eines Schichtstapels, wie Schichtdicken, 
Brechungsindizes oder Absorptionskoeffizienten, durch Vergteich eines 
von einem Ort des Schichtstapels aufgenommenen optischen 
MefJspektrums mit einem anhand vorgegebener optischer 
Parameterwerte berechneten Analysespektrum und Optimierung des 
berechneten Analysespektrums auf das Meftspektrum hin. Hierbei wird 
vorgeschlagen, das aufgenommene Melispektrum anhand von 
Kurvenformparametern, die das MelJspektrum charakterisieren und aus 
diesem ermittelt werden, zu klassifizieren und diese Kurvenformparameter 
mit entsprechenden fur bekannte Schichtstapel berechnete 
Kurvenformparametern von Spektren zu vergleichen, urn (Ausgangs-) 
Werte oder Wertebereiche fur die zu bestimmenden optischen Parameter 
zu ermitteln, anhand derer das oder die Analysespektren zum Vergleich 
mit dem MefSspektrum errechnet werden. Die Erfindung erlaubt eine 
drastische Reduzierung von Rechenkapazitat und -zeit 
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